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We develop a generalized version of heavy-baryon chiral perturbation theory to describe pion- 
nucleon scattering in a kinematic domain that extends continuously from threshold to the delta- 
isobar peak. The P-wave phase shifts are used to illustrate this framework. 
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1. Introduction 

A prominent feature of pion-nucleon {nN) scattering is the delta resonance, A(1232), a peak 
in the elastic cross section at the center-of-mass (CM) energy wia = mj^ + 5 ~ 1230 MeV, where 
8 ~ 290 MeV is the nucleon-delta mass splitting A resonance can be studied by considering the 
unitarity and analyticity of the 5 matrix; however, the accuracy is hard to improve systematically 
with these general principles alone. Our goal here is to investigate TlN scattering from threshold up 
to the delta resonance in an effective field theory (EFT) (for more details, see Ref. iQ]). 

Following several seminal papers [^, EFTs have been developed as model-independent ap- 
proximations to low-energy strong interactions, which can be systematically improved by a series 
in powers of Q/Mqcu, where Q refers generically to small external momenta and Mqcd 1 GeV 
is the characteristic QCD scale. For reviews, see, for example, Refs. |Q §]. Chiral perturbation 
theory (ChPT) specializes in processes involving at most one nucleon [Q]. ChPT with only pion 
and nucleon fields has been extensively applied to near-threshold nN scattering resulting in a 
perturbative expansion in powers of 2/5 and mji/5, which converges slowly a.smji/5 ~ 1/2. The 
convergence can be improved with an explicit delta field. The explicit delta in nN scattering within 
standard ChPT has been explored and demonstrated in a fully consistent calculation [||]. 

Nevertheless, the perturbative nature of standard ChPT makes it impossible to describe the 
delta resonance, a non-perturbative phenomenon. A non-perturbative treatment of the delta within 
ChPT was considered in Ref. ; however, a systematic resummation did not exist until the seminal 



work of Ref. [10], where it was justified by a power counting based on three separate scales mj^ <^ 



5 <C Mqcd, and this idea has been applied to various electromagnetic reactions in the delta region. 



but for nN scattering few results have been published [|ll|]. 

We employ a power counting developed for generic narrow resonances [12], in which there are 
only two scales Mio ^ 5 ^niji and Mhi ~ Mqcd- Thus the EFT expansion of the nN scattering am- 
plitude pursued here is in powers of Q/M\^ and Mio/My. The kinematic region under consideration 
spans over both threshold and the resonance. 



2. Effective Lagrangian 

To establish the notation, we review how the effective Lagrangian is constructed (for more 
details of building chiral Lagrangian, see e.g., Refs. [|, |, 0, [T|]). The effective La grangian 
should inherit the symmetries of QCD: Lorentz invariance, (approximate) two-flavor chiral sym- 
metry {SU {2)i X SU (2)s), parity, time-reversal invariance, and baryon-number conservation. 

In the kinematic region where the EFT holds, external momenta are much smaller than the 
nucleon mass, Q <^ niN, and thus Lorentz invariance can be fulfilled perturbatively in powers of 
Q/rriN. One can start with a relativistic Lagrangian using the Rarita-Schwinger field for the delta, 
and then reduce from it its nonrelativistic version 11]. This way, however, extra effort needs to 



be taken in order to control the spurious spin- 1/2 sectors of the Rarita-Schwinger field. We employ 
another approach that starts with heavy-baryon fields N for the nucleon and A for the delta, which 
are, respectively, two- and four-component spinors in spin and isospin spaces. Eventually, the 
effective Lagrangian only has the baryon degrees of freedom that represent forward propagation. 
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The crucial ingredient in this approach is to develop an order-by-order Lorentz transformation, by 



which one can constrain the coefficients of the rotation-invariant operators QlSj |17| ]. 

Due to the presence of the delta field A, one needs 2x4 matrices S in spin space to make a 
three-vector bilinear, and Q.ij a three-tensor. Similar transition matrices, T and "Eab, can be 
defined in isospace. 

The chiral-invariant operators are isoscalars that are made of pion covariant-derivative = 
D^^d^ll/lfji with D = 1 + n^/4f^, N, A, and their covariant derivatives, for example, i^^A = 
d^+t('^) -E^) A with Ef, = i7t/f^xD^. 

We use the so-called chiral index v |^ to organize the operators of the effective Lagrangian 
v = d + m + n§+ f /2 — 2, where d, m, ng, and / are the numbers of derivatives, powers of nijc, 
powers of 5, and fermion fields, respectively. In constructing the Lagrangian, we use integration 
by parts and field redefinitions to remove time derivatives on baryon fields except for the kinetic 



terms. The Lagrangian terms with the two lowest indices are given by [15] 



= 2f^D^-^ml7c'^+N'^i%N + gAN'^tdN--b 

+A^ {i% - 5) A + 4g^A^?(2)5(i)A. -D + ZiA (n^TSA + H.c^ ■ b + 



and 



2mN 

while the next-higher index yields 



^ ^N^^^N + A^^^a] - — iiN^TS-&A + H.c.]-Do + - 



(2.1) 



(2.2) 



^(2) 



5 + -?o hi 



+ 



2m2 
8m?, 



2m2 



(n'^TS^^A-N'^TS-^^A^ +H.C. 



■ D 



dimN^TS- &A + 3N^TSi&m^ + 2£ijiN^Tai^Si;A\ +H 



.c. 



+di (N'^TSA + H.c)j--^^b + d2^ {N^TSA + Hx.y-b + --- (2.3) 



Here, gA (Sa) is the v = axial-vector coupling of the nucleon (delta) and Ha (t/i,2) is (are) the 
V = (v = 2) ;rA^A coupUng(s). 



3. Power counting 

When the CM (heavy-baryon) energies E are much below the delta peak, the power counting is 
standard [||, ^, |]] with the simple generalization that 5 counts as Q. The contribution of a diagram 
with A nucleons (here A = 1), L loops, and Vj vertices with chiral index v, is proportional to QP, 
with 

p = 2-A + 2L + £v,V;. (3.1) 

i 

However, in the small region spanning the delta peak whose size is of the leading-order (LO) 
delta self-energy, — 5| ~ Z^"' = ^^{Q^ /Mq^^), a resummation is needed in one-A-reducible 
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Figure 1: Contributions to nN scattering up to order Q : (A) Q ' pole diagram; (B) pole dia grams; (C) 
pole diagrams; (D)&(E) tree dia grams, of which (E) apply to both regions. E^"' is the n-th order delta 
self-energy, V^j"' the nNA vertex function, and zjj"' (Z^"') the nucleon (delta) field renormalization constant. 



diagrams because one insertion of £^°' and the bare delta propagator contributes 
5) = 1^(1). The resummation thus amounts to a dressed propagator 



(1): Zf/[E 



sf{E) = E-5 + I>:'{5) 



(0), 



-1 



(3.2) 



which scales as Mq^^q/Q^. This is an enhancement of two powers over the generic situation. As a 
consequence, in one-A-irreducible diagrams the standard ChPT power counting (3J) still applies; 
dressed delta propagators only need to be included in one-A-reducible diagrams. We thus arrive 
at a new power counting for one-A-reducible diagrams within a narrow window around the delta 
peak, 

p =2-A-ln^ + 2L + YViVi , (3.3) 

i 

where «a is the number of dressed delta propagators. This is the non-electromagnetic version of p 
derived in a slightly different power counting in Ref. [10]. Diagrams up to are listed in FIG. |l[ 
It seems that the two different power-counting schemes, which are applicable in two different 
regions, would lead to an EFT amplitude in the form of a piecewise function in the energy. Even 
worse, separating these two regions is somewhat arbitrary. A piecewise EFT is actually unnecessary 
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if we enforce the pole diagrams even in the off-the-pole region, which is equivalent to shifting 
a subset of higher diagrams into lower orders, i.e., a rearrangement of diagrams. This sort of 
rearrangement still retains the essence of the original power counting as long as one does not claim 
a higher accuracy by doing so. 



4. ;rA^-scattering T matrix 



The partial-wave T matrix is related to the phase shifts, in the channel with total angular 
momentum j, orbital angular momentum /, and isospin t, by 



Tju{E) = -i{tx^[2ieju{E)\-\] . 



(4.1) 



In the following we will use a more conventional notation for a specific partial wave: ht^ij- For 
example, P\t, refers to the / = 1 (P wave), f = 1/2, and j = Ji/l. 

Here the exact relation between E and the CM momentum k, E = {mj^ + k^Y^'^ + (m^ + 
mliY^^ —mj^, is assumed, meaning that certain trivial, kinematic k/mp^ terms are resummed — 
what we refer to as semi-resummation. A strict heavy-baryon expansion can be readily obtained 
afterwards. 

At LO (2^^) there is only a pole diagram, FIG. |l|(A), which contributes only to the Pj,^ wave, 

y(o)(5) 



7.L0 



£-5 + /7(o)(5)/2 



1 + ^ 



M, 



Q_ 

QCD 



(4.2) 



where 
7(0) (5) 



h\ /c-2 2\7 ri c- / /c-2 2n //^ sill '^ + ^/^N+{8'^ -m\)/2ml. 



24nf^ 

The NLO (2") amplitude has the same form as LO, 

y(0)(5) + ^(i)(5) 



rNLO 

'P33 



£-5 + / [7(0) (5) + 7(0(5)] /2 



1 + - 



(4.3) 



(4.4) 



However, 7^0(5) vanishes in the CM frame when we do not expand kinematic relations in powers 
of 

Summing up the pole (FIG. |l|(C)) and tree (FIG. [l](E)) diagrams, one first finds the NNLO 
amplitude in the P33 channel. 



T^NNLO 



r(£) 



E-5 + ir{E)/2 



[l + iTB{E)] + TB{E) + ff{Ti 



LO 
^■33 



M, 



(4.5) 



QCD^ 



where 



Tb{E) 



(§1 , 1 hi 
671/2 36E + d 



and r(E) 



[m 



E + mj^ 



with 



k^ 



hA 



mt 



(4.6) 



(4.7) 
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Figure 2: phase shifts (in degrees) as a function of Wcm (in MeV), the CM energy including the nucleon 
mass. The EFT strict heavy-baryon expansion at LO (NNLO) is represented by the red dashed (black solid) 
line. The NLO curve coincides with LO. The light-blue band outlines the estimated systematic error of the 
NNLO curve. The green dots are the results of the GW phase-shift analysis [|l3]]. Points marked by a red 
star (black square) are inputs for LO (NNLO). 



where 



+ 



i-/il(l3+15x2)+^gf 



In 



Vi+x2 + i 



(4.8) 



and kg satisfies 5 = (m^ + kg)^^^ + (m| + k^y/'^ — nif^. Other channels are easy to calculate from 
the one-A-irreducible tree diagrams in FIG. |I](E). For the remaining P-wave channels, 



rNNLO 
'A3 



T^NNLO 



— T 



NNLO 



1 hi 



nnfl\E 9E + 5 



1 + 1 



Q 



M, 



QCD 



(4.9) 



5. P-wave phase shifts 

A number of low-energy constants (LECs) can be determined from other processes, such as 
pion decay and neutron decay. We adopt the following values: niji = 139 MeV, nif^ = 939 MeV, 
gA = 1 -26, and fjc = 92.4 MeV. Our strategy of fitting is to determine the free parameters, 5, Ha, and 
X from the P33 phase shifts around the delta peak and then predict the phase shifts at lower energies 
in all P waves. Shown in FIGs. | and | are the EFT curves (strict heavy-baryon expansion used) 
fitted to the partial-wave analysis (PSA) by the George Washington (GW) group [[T3|]. The PSA 
points used to determine the free parameters are explicitly marked. Based on the power counting we 
use, systematic errors of the EFT curves can be estimated, shown in FIGs. ^ and ^ as light-colored 
bands. 

The LECs extracted from the P33 fit are given in TABLE [ij One can estimate the errors in the 
NNLO values as the variation in each LEC within which the NNLO P33 curve in FIG. ^ roughly 
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Figure 3: Predicted phase shifts (in degrees) in the P\^, P-n, and P\i channels as functions of Wcm (in 
MeV), the CM energy including the nucleon mass. LO and NLO vanish in these channels; NNLO EFT 
results in the strict heavy-baryon expansion are given by the black solid lines. The light-blue bands outline 
the estimated systematic errors of the NNLO curves. The green dots are the results of the GW phase-shift 
analysis [jl^. 



Table 1: Low-energy constants extracted at LO, NLO, and NNLO from the fits using the strict heavy-baryon 
expansion. 



8 (MeV) 
LO NLO NNLO 


hA 

LO NLO NNLO 


NNLO 


293 293 321 


1.98 4.21 2.85 


0.046 



stays within the eiTor band. This way we find 5/MeV, Ha, and x to be within ~ ±4, ±0.30, and 
±0.030, respectively, of the NNLO values in TABLE [I[ 



6. Summary 

We have extended standard ChPT to deal with the non-perturbative delta resonance in an EFT 
framework. The delta is treated as a nonrelativistic particle from the beginning, rather than being 
represented by the Rarita-Schwinger field. 

Like other EFTs that deal with non-perturbative phenomena, ours captures the non-perturbative 
structure in LO. Subsequently, the power counting leads to a systematic, perturbative improvement 
beyond LO. We applied this power counting to low-energy nN scattering, where we built the am- 
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plitudes up to NNLO. We fitted our P-wave amplitudes to the phase shifts given by Ref. [13]. With 
just three free parameters, we obtained a good fit in the ^33 channel. 

The EFT approach presented here also provides the basis for a model-independent, unified 
description, from threshold to past the delta resonance without discontinuity, of reactions involving 
other probes and targets, including nuclei. 
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